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Dilution effects on the long-range ordered state of the doubly degenerate eg orbital are investigated. Quenched 
impurities without the orbital degree of freedom are introduced in the orbital model where the long-range order 
is realized by the order-from-disorder mechanism. It is shown by the Monte-Carlo simulation and the cluster- 
expansion method that a decrease in the orbital-ordering temperature by dilution is remarkable in comparison 
with that in the randomly diluted spin models. Tilting of orbital pseudo-spins around impurity is the essence 
of this dilution effects. The present theory provides a new view point for the recent resonant x-ray scattering 
experiments in KCui_ i.ZnvF3. 

PACS numbers: 75.30.-m, 71.23.-k, 71.10.-w, 78.70.-g 



Impurity effects on the long-range ordered state are one of 
the attractive themes in recent study of correlated electron sys- 
tems 0]. A small amount of non-magnetic impurities dramati- 
cally destroy the superconductivity in cuprates, and induce the 
antiferromagnetic long-range order in the low-dimensional 
quantum spin-liquids. Doped impurities also cause striking 
effects on the charge- and orbital-orders; substitution of Cr 
ions for Mn in colossal magnetoresistive (CMR) manganites 
immediately destroys the charge-orbital orders Origin 
of CMR itself is studied from the view point of randomness 
and/or percolation [3J. 

Recently, the dilution effects in KCUF3 by substituting Zn 
for Cu is reported experimentally by Tatami etal. 01 • A 
Cu^+ ion in the cubic-crystalline field shows the flij^g con- 
figuration which has the eg orbital degree of freedom. The 
Cu ions in the perovskite crystal form the three-dimensional 
(3D) simple-cubic (SC) lattice, and exhibit the long-range or- 
bital order (OO) in room temperatures, where the d^i_,i- and 
d^2_^2-\\ke orbitals are aligned with momentum Q = (tt, tt, tt). 
Since the five d orbitals are fully occupied in Zn^+, substitu- 
tion of Zn for Cu corresponds to dilution of orbital. The res- 
onant x-ray scattering (RXS) studies in KCui_^Znj^F3 reveal 
that a decrease in the orbital-ordering temperature (Too) by 
dilution is remarkable in comparison with the randomly di- 
luted magnets, and OO disappears around x — 0.5, as shown 
in the inset of Fig. 1. These observations may not be ex- 
plained by the conventional percolation scenario; the site- 
percolation threshold in 3D SC lattice is Xc{= 1 — Pc) =0.69 
which is applicable well to the several diluted magnets such 
as KMni_,Mg,F3ill3. 

We examine, in this Letter, the dilution effects on the long- 
range order of the eg orbital degree of freedom. As well 
known, the doubly degenerate eg orbital is treated by the 
pseudo-spin (PS) operator with a magnitude of 1/2; T, = 
\Y,yyi ,sd]y^Gyyidiyis whcrc t/jyj Indlcatcs the annihilation op- 
erator of a hole with spin s and orbital 7 at site /, and a is the 
Pauli matrices. A shape of the electronic orbital is described 
by an angle of the PS vector as \Q) — cos{Q /2)\d-^,2_yi) + 
5i\\{d /2)\d^2_yi). For example, 0=0, 2lt/2> and 4;r/3 cor- 
respond to the d^^2_y2, d^^2_y2, and d^^2_y2 orbitals, respec- 
tively, and 9 = n/3,n and — 7r/3 to d^2_^2, d^2_y2, and d^2_^2. 
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FIG. 1: Impurity concentration dependence of the normalized 
orbital-ordering temperature Too (^)/ Too ~ 0). The calculated 
results by the MC method are shown by the closed- and open-red 
squares (see text), and those by the CE method for the one- and two- 
site clusters are shown by the bold- and broken-blue lines, respec- 
tively. The results by the CE method where the PS operator is treated 
as a classical vector are shown by the one-point chain line (blue). 
For comparison, Tj<j in the XY (reverse-triangles) and Ising (stars) 
models by the MC method, and that in the Heisenberg model (dot- 
ted line) by the CE method are also shown. The inset shows the Zn 
concentration dependence of Too in KCui_j:Znj:F3 obtained by the 
RXS experiments ||4(|. 



respectively. As shown below, the Hamiltonian for the diluted 
orbital system is described by the PS operators and seems 
to be on the same footing with the spin models. However, 
the obtained results are qualitatively different from the di- 
luted magnets. The main results are shown in Fig. ^ The 
impurity-concentration x dependences of Too, calculated by 
the two different methods, show rapid decreases in compar- 
ison with that in the spin models. The results can explain 
the experimentally discovered anomalous dilution effects in 
KCui_;,ZnvF3 H. 

The model Hamiltonian adopted here describes the orbital 
interaction between the nearest-neighbor (NN) Cu sites in a 
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3D SC lattice: 

^^-27^T/Tje,e;, (1) 

where J{> 0) is the interaction, and (ly) indicates a pair of 
the NN sites along the direction l{=x,y,z). The operator t/ 
depends explicitly on the bond direction and is defined by a 
linear combination of the PS operators as 

t:' = cos {-J-'^i^ Ti- + sin (^^"'^ ^-f (2) 

with (nj-,nj,,n,) = (1,2,3). The quenched impurities without 
the orbital degree of freedom are introduced by e,- taking one 
and zero for Cu and Zn, respectively. 

The Hamiltonian in Eq. ([0 without impurities (e, — 1 for 
VO does not concern an origin of the interaction, i.e. the elec- 
tronic and/or phononic interactions. In the electronic view 
point, this is derived by the generalized-Hubbard model with 
the eg orbital degree of freedom. Through the perturbational 
expansion with respect to the NN electron transfer f, the 
spin-orbital superexchange model is obtained 1^ ; = 
-2/1 1(,7)(i + S, • S;)(| - T/rj) - 2y2l(,;)(i - S,- • Sj)il + 
r/rj + t' + fj). Here, S; is the spin operator at site / with 
magnitude of 1 /2, and Ji (= t^/{U - 3/)) and J2{= t^/U) are 
the superexchange interactions with the on-site intra-orbital 
Coulomb interaction U and the exchange interaction /. Since 
the Neel temperature (7n) in KCui-^ZniFg is far below Too in 
a whole range of x 1 4|, taking S, • = is a good assumption, 
and Eq. Q with / = |7i — is obtained. In the phononic 
view point, the orbital model is derived based on the coopera- 
tive Jahn-Teller (JT) effects. Start from the linear JT coupling 
Hamiltonian ,yfjT ~ gY^i,m={x.7)QimTim with the vibrational 
modes 2,„, [m = x,z) in a Eg octahedron, and the lattice po- 
tential K between NN Cu-E bonds. After rewriting 2»n by the 
phonon coordinates qt, we obtain Eq. Q with J — 2g^/ {9K) 
by integrating out qk A sign of J is positive in both 

the two processes. A unique aspect of the orbital model to 
be noticed here is that the explicit form of the interaction de- 
pends on the bond direction /; the interactions are Tlv7}j. for 

I = z, and [- i T„ + (- ) ^ Ti,] [- \ Tj, + (- ) ^ Tj,] for l^x{y). 
When we focus on one direction I, the staggered alignment of 
the d^i2_^2 and t/„,2_„2 orbitals with {l,m,n) = (jc,y,z), {y.,z.,x) 
and {z,x,y) is favored. In a 3D SC lattice with and without 
impurities, the stable orbital configurations are non-trivial. 

We have attacked this issue by the numerical approach to- 
gether with the analytical one, i.e. the classical Monte-Carlo 
(MC) simulation, and the cluster expansion (CE) method. In 
the MC method, the PS operator is treated as a classical vector 
defined in the Tx — plane. The MC calculations have been 
performed for the cubic Lx Lx L lattice (L = 10 ^ 18) with 
the periodic-boundary condition. Eor each sample, 30,000MC 
steps are spent for measurement after 8,000MC steps for ther- 
malization. The physical quantities are averaged over 20 ^ 
80 MC samples at each parameter set. We adopt the CE 
method proposed in Ref. il 111 , where the fluctuations of the 



effective fields are determined with the order parameter, self- 
consistently. Even in the two-site clusters, this CE method 
provides good values for Xc and the Curie temperature 7c at 
x = for Ising and Heisenberg models 11 ill . 

First, we show the results without impurities. It is known 
that, in the orbital model [Eq. Q] at x = 0, there is a 
macroscopic number of degeneracy in the mean-field (ME) 
ground state 

O El El Eld These are classified into 
the two; (i) One of the MF solutions is the staggered-type 
OO with two sublattices, termed A and B, and the momen- 
tum Q = {n,n, 7z). The orbital angles in the sublattices are 
[Sa/Ob) = (0/^^ + ^) for any value of 9. Such continuous- 
rotational symmetry is unexpected from Eq. Q. (ii) Con- 
sider an OO with Q = {n,n,n) and (0a/0b) = {d(,/d(, + K), 
and focus on one direction in a 3D SC lattice, e.g. the 
z direction. The MF energy is not changed by changing 
[Qq/Qq + ti) {—9o/~9Q~n) in each layer in the xy plane. 
Both the types of degeneracy are understood from the mo- 
mentum representation of Eq. ([Q: — Y.k Wk-^i^)Wk with 
i/k — [Tiki^xk] and the 2x2 matrix 7(k). By diagonahzing 
7(k), we obtain the eigen values /±(k) = 27[— c,- — Cy — ± 

^c^ + Cy + c-^ — CxCy — CyC. — Cj-c-i] withc/ =cos(;r^//a) JitIi . 

y_(k) has its minimum along (;r,;r,;r) — {Q,n,7z) and other 
three-equivalent directions. A liftin g of the degeneracy in the 
orbital model is discussed in Refs. 13 Ell by the spin-wave 
analyses. Here we show the lifting of the degeneracy by the 
MC method. In Fig.|2l(a), number of data obtained in the sim- 
ulation is plotted as a function of the orbital angle in the 
staggered-type OO with {9a/ Ob) = {9 /O + 7z). It is shown 
that the type-(i) orbital degeneracy is lifted; 9 is confined to 
the three discontinuous values of 0, 2;r/3 and 4n/3. The inset 
shows number of data as function of the orbital-order param- 
eter Mo at Q = (n, n, %) defined as = L«=.T,z(7Qm7'-Qm) 
with rk,„ ](?L, e''""' 7;,„. Most of the data show Mq = 0.5, 
implying that the ordering pattern is of a staggered type, and 
the type-(ii) degeneracy is also lifted. 

We have performed the finite-size scaling analyses in the 
MC simulation to determine Too. The correlation length t, is 
calculated by the second-moment method in the PS correla- 
tion function for several sizes L. In Fig.|2jb), we demonstrate 
the scaling plot for ^/L as a function of {T - Too)L}I'' at 
X = 0. The scaling analyses work quite well for L =10, 12 
and 14. Too and the exponent v are determined through the 
least-square fitting by the polynomial expansion and obtained 
as Too /J = 0.344 ± 0.002 and v = 0.69 - 0.8 1, although sta- 
tistical errors are not enough to obtain the precise value of 
V. Even in the diluted case [see the inset of Fig. |2jb)], the 
precision is enough to determine Too[= (0.248 ±0.003)7 for 
X ~ 0.15 ]. Beyondx = 0.15, the scaling analyses do not work 
well. 

In Fig.|3ja), the temperature dependence of the normalized- 
order parameter -|J-^Mo[= oto(x, T)] are presented. First, fo- 
cus on the region of x < 0.15. As expected, mo{x = 0,T) 
abruptly increases at Tqo and is saturated to 0.5 at T ^ 0. By 
doping impurity, mo(x, T) does not reach 0.5 even far below 



3 




T/J 



FIG. 2: (a) A MC data distribution for the orbital angle 9 in the 
staggered-type orbital order at x = 0. The inset shows a data distri- 
bution of the order parameter Mq at x = 0. (b) The finite-size scaling 

for the correlation length £, /L as functions of L'/*' (T - Too) /J at FIG. 3: (a) Temperature dependence of the normalized order param- 
x = 0. The inset shows the scaling plot at x = 0.15. eter mo(x, r)[= Mo/(l -x)] for various x and L. (b) A snapshot in 

the MC simulation for the PS configuration at x = 0.1. Open circles 
indicate the impurities, (c) A snapshot at x = 0.3. 



Too, and mo{x 7^ 0, T ^ 0) gradually decreases with increas- 
ing X. Although the system sizes are not sufficient to estimate 
mo{x, T) in the thermodynamic limit, mo{x ^ 0,T ^ 0) does 
not show the smooth convergence to 0.5 in contrast to the di- 
luted spin models. Beyond x = 0.15, the situation is changed 
qualitatively; in spite of the fact that mo{x, T) grows around a 
certain temperature (e.g. ~ 0.187 at x = 0.2), mo{x, T ^ 0) 
becomes small abruptly and decreases weakly with increas- 
ing L. Anomalies in the specific heat and the susceptibility 
around this temperature become weak, and some of the phys- 
ical quantities, e.g. A^o = Lm=x,z(7Qm)^ with Q = {n,n,n), 
depend on the initial states in the MC simulation. It seems 
that the orbital correlation grows up below this temperature, 
but the long-range order does not. We suppose that a possi- 
ble orbital state in this region is a short-range ordered state 
or a glass state 1 18], although we do not identify this phase 
coiTectly at the present stage. Thus, in jc > 0.15, a tempera- 
ture where dmo{x, T)/dT in L = 18 takes a maximum is in- 
terpreted to be the cross-over and/or glass-transition tempera- 
ture. Above X = 0.3, estimation of Too becomes severe, be- 
cause of the weak temperature dependence of mo{x, T). Thus, 
as a supplementary information of the ordering temperature, 
we calculate a temperature Too where dNo/dT takes a maxi- 
mum, and consider its full width at half maximum as an error 
The initial state in the calculation for Tqq is assumed to be the 
ordered state, although the initial-state dependence of Too is 
much weaker than that of the A^o amplitude. 

The X dependences of Too (and Too) presented in Fig.fjare 
obtained by the MC (filled-red squares) and CE (blue lines) 
methods introduced above. For comparison, the x depen- 
dences of 7n in the spin models are plotted. Tbos are also 
shown by the open-red squares. Starting to dope impurities 
in the orbital model, decrease of Too is more remarkable than 



that of 7n in the spin models. Around x = 0.15, the clear 
transition to the long-range order becomes obscure, explained 
above. Toos obtained by A^o are close to Too below x = 0.2, 
and decrease smoothly up to x = 0.3, although the critical x, 
where Too disappears, is not determined due to large statis- 
tical eiTors. A rapid decrease of Too is also obtained by the 
CE method. It is seen that difference between the calculated 
Too — X curves for the cluster sizes At = 1 (bold line) and 2 
(broken like) is within a few percent. Too in the CE method 
disappears around x = 0.5 being much smaller than Xc = 0.69 
for the 3D SC lattice. We also calculate Too in the CE method 
where the quantum PS operators are replaced by the classical 
vectors. The calculated results plotted by the one-point chain 
line are close qualitatively to the Too (Too) curve obtained 
by the classical MC method. 

Now we introduce the physical picture of the diluted 
orbital-ordered state. The real-space configurations of the 
orbital PSs help us to understand the calculated Too ~ x 
curve. The MC snapshots of the PSs on a plane are shown 
in Figs. |3jb) and (c) for x=0.1 and 0.3, respectively. The 
staggered-type OO with the orbital angles (Oa/Ob) = (O/tz) 
is seen in the back-ground of Fig. |3lb). At the neighboring 
sites of the impurities indicated by the open circles, tiltings of 
the PS vectors are observed. Disturbing the PS configuration 
from (Q/tz) becomes violent at x = 0.3. These observations 
are caused by the local-symmetry breaking by dilution. Con- 
sider the orbital state at a neighboring site of an impurity on 
the X axis. On account of the impurity, one of the interactions 
along X vanishes. Since the interaction depends on the bond 
direction explicitly in the orbital model, the PS at this site tilts 
to gain the interaction energies for other five bonds. This is 
the essence of the diluted orbital systems and is in contrast to 
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FIG. 4: Impurity concentration dependence of the normalized 
orbital-ordering temperature Too{x) /Too{x = 0) with the higher- 
order JT coupling. The squares, circles and triangles indicate 
Too{x)/Too{x = 0) at j' /J =0, 0.3 and 3, respectively. Toos in- 
dicated by the open and filled symbols are defined in the text. For 
comparisons, of the XY model (reverse triangles), and Tqq at 
j' /J=0 by the CE method (bold line) are also plotted. 



the conventional diluted spin models. To check roles of this 
PS tilting on the anomalous Too —x curve, we have performed 
the MC calculations in the model where the available orbital 
angles in the T^ — T^ plane are restricted to a few states. The 
obtained Too — x curve is similar to that in the spin models; 
the tilting of PSs is confirmed to be the origin of the rapid 
decrease of Too. 

Finally, to compare the present theory with the experiments 
in KCui_;rZn^F3, more directly, we introduce the higher-order 
JT coupling. From the RXS experiments |4] and the lat- 
tice distortion in KCuFa Il9l] . the OO in KCuFa is expected 
to be the dy2_,2/d^2_-^2-type; the cant-type OO with the or- 
bital angles (07- 0) (0 ^ 7t/3) and Q = To re- 
produce this type of OO, the higher-order JT coupling is 
required, J^yr = g"L, { {Q\ - QI) Ti, - IQi.Qi.Ti,) which 
provides the anisotropy in the PS space. Here we treat ap- 
proximately this term based on the theory of the cooperative 
JT effects. By rewriting by the phonon coordinates 
Qk, and integrating out qk, the interaction between orbitals 
at the different three-neighboring sites is obtained: — 
J' I(/^> { {TizTjz + Ti,Tj,) Tfe - ITiJjJk, } £,ejek- Here, / = 
{I6g' g^) / {9K'^) is the coupling constant and {ijk) indicates 
a sum of the neighboring three sites. The MF energy of 
is proportional to =j-cos30 implying the anisotropy. A value 
of g' for a Cu^+ ion was estimated by the adiabatic potential 
baiTier in a molecule |20| and corresponds to J' being about 
0.37. We have performed the MC calculations in the model of 
+ and observed that the cant-type OO at x = is real- 
ized. The doping dependences of Too (and Too) with includ- 
ing are presented in Fig.|3for7'/7 — 0.3 and 3. The filled 
and open symbols are for Too and Too, respectively, which are 



obtained by the same ways with those in Fig. 1 . With increas- 
ing /, the 7bo(7bo) -~x curves approach to those for the spin 
models, because the anisotropy suppresses the tilting of PSs. 
However, we confirm that the rapid reduction of Too by dilu- 
tion survives even in the calculation with the realistic value of 
/, and is consistent with the RXS experimental results shown 
in the inset of Fig. 1. 

In summary, we have investigated the dilution effects on 
the long-range order of the eg orbital degree of freedom. We 
confirm, by both the MC and CE methods, that Too decreases 
rapidly with doping, in comparison with the diluted magnets. 
The tilting of the orbital PSs around impurities, which is dis- 
tinguished qualitatively from the spin models, is the essence 
for the rapid decrease of Too . The present theory provides 
a new view point for the recent experiments in a diluted or- 
bital systems of KCuj _;fZn^F3 . A broad peak profile observed 
by the RXS experiments ||4|] may be attributed to the orbital 
tilting around impurities. Observations of the tilting by the 
scanning-tunneling microscope and/or the scanning-electron 
microscope are a direct check for the present results. 
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